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Abstract 

We discuss the response of a Bose-Einstein condensate to a change in the scattering length, which 
is experimentally realized by tuning the magnetic field near a Feshbach resonance. In particular, we 
consider the collapse of the condensate induced by a sudden change in the scattering length from 
a large positive to a small negative value. We also consider the condensate dynamics that results 
from a single pulse in the magnetic field, due to which the scattering length is rapidly increased 
from zero to a large value and then after some time rapidly decreased again to its initial value. We 
focus primarily on the consequences of the quantum evaporation process on the dynamics of the 
Bose-Einstein condensate, but also discuss the effects of atom-molecule coherence. 
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I. INTRODUCTION 



From single-channel scattering theory it is well known that the collisional cross section 
changes dramatically if the energy of the incoming particles is close to the energy of a long- 
lived bound state in the interaction potential. In particular, the magnitude of the s-wave 
scattering length a of an attractive potential well diverges as the depth of the potential well 
is increased such that a new bound state enters the well A similar behavior occurs in the 
case of a Feshbach resonance, when the energy of the two particles in the incoming channel is 
close to the energy of a bound state in a closed channel 0. In the case of collisions between 
alkali atoms, the coupling between the two channels is provided by the hyperfme interaction. 
Due to the spin flips involved in this interaction, the difference in energy between the bound 
state and the continuum, the so-called detuning, is adjustable by means of a magnetic bias 
field. Feshbach resonances were first predicted theoretically [|, but have now also been 
observed experimentally, in various atomic species [H, H, H, B]- As a result the experiments 
with magnetically-trapped ultracold atomic gases, where the s-wave scattering length fully 
determines the interaction effects, have an unprecedented high level of control over the 
interatomic interactions. In this paper we focus on 85 Rb in the / = 2,m/ = —2 hyperfme 
state, which has a Feshbach resonance at a magnetic field of Bo « 154.9 (G)auss ||. Near 
the resonance the scattering length, as a function of magnetic field, is given by 

a{B) = a h Jl - ——) . (1) 
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The resonance is characterized by the width AB « 11.6 G and the off-resonant background 
scattering length a^ g ~ — 450a , with a the Bohr radius. In Fig. [I] the scattering length is 
shown as a function of the magnetic field. Clearly, it can be adjusted experimentally from 
large negative values, to large positive ones. Moreover, at a magnetic field of B « 166.5 G 
the scattering length vanishes, and the gas behaves effectively as an ideal Bose gas. 

With this experimental degree of freedom it is possible to study very interesting new 
regimes in the many-body physics of ultracold atomic gases. The first experimental appli- 
cation was the detailed study of the collapse of a condensate with attractive interactions. 
In general a collapse occurs when the attractive interactions overcome the kinetic energy of 
the condensate atoms in the trap. Since the typical interaction energy is proportional to 
the density, there is a certain maximum number of atoms above which the condensate is 
unstable |10|, [□], [T^, [13], |i~4"l . In the first observations of the condensate collapse by Bradley 



et al. fi~5H , a condensate of doubly spin-polarized 7 Li atoms was used. These atoms have 
a magnetic-field independent scattering length of a ~ — 27ao- For the experimental trap 
parameters this leads to a maximum number of condensate atoms that was so low that non- 
destructive imaging of the condensate was impossible. Moreover, thermal fluctuations due 
to a large thermal component make the initiation of the collapse a stochastic process, thus 
preventing also a series of destructive measurements of a single collapse event [16|, [L7], |i~8| . 
A statistical analysis has nevertheless resulted in important information about the collapse 



process. Very recently, it was even possible to overcome these problems [19 . 

In addition to the experiment with 7 Li, experiments with 85 Rb have been carried out. 
In particular, Roberts et al. fl2(J also studied the stability criterion for the condensate, and 



Donley et al. jn[] studied the dynamics of a single collapse event in great detail. Both 
of these experiments make use of the above-mentioned Feshbach resonance to achieve a 
well-defined initial condition for each destructive measurement. It turns out that during 
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a collapse a significant fraction of atoms is expelled from the condensate. Moreover, one 
observes a burst of hot atoms with an energy of about 150 nK. Several mean- field analyses of 
the collapse, which model the atom loss phenomenologically by a three-body recombination 
rate constant, have offered some theoretical insight ^3|, |25|, |26], fT7| . However, the 
physical mechanism responsible for the explosion of atoms out of the condensate and the 
formation of the noncondensed component is still largely ununderstood at present. 

A second experimental application of the Feshbach resonance has been implemented by 
Claussen et al. [J2H]. Starting from the noninteracting limit, the scattering length is made 
to jump very fast back and forth to a large and positive value. Surprisingly, in this case one 
also observes loss of atoms from the condensate, as well as a burst of hot atoms. Even more 
surprising is the fact that the amount of atoms expelled from the Bose-Einstein condensate 
decreases with time during some intervals, excluding a theoretical explanation in terms of a 
loss process characterized by a rate constant. 

As a third application Donley et al. |3l| have conducted an experiment where two trape- 
zoidal pulses in the magnetic field were applied. As a function of the time between the 
two pulses an oscillatory behavior in the number of condensate atoms is observed, which 



is attributed to coherent Rabi oscillations between atoms and molecules \E2, 33, EM. In 



this paper we will not consider this experiment, but instead focus on the first and second 
experimental applications. In particular the single-pulse experiment has not received much 
attention theoretically, even though an understanding of this experiment seems an essential 
first step in the theoretical treatment of the recent two-pulse experiment. Therefore, the 
discussion of this most recent experiment will be postponed to a future publication. 

In a previous paper we have considered the loss of atoms by means of elastic two-body 
collisions, in the situation where the condensate collapses |2I5j . However, the mechanism put 
forward in this paper is much more general. In particular, the mechanism should also be 
relevant for the above mentioned single-pulse experiments. The main goal of this paper is 
to present the theory behind it in great detail. Since the mechanism is able to describe loss 
from a condensate at zero temperature we will hereafter refer to it as quantum evaporation. 
The two-pulse experiments of Donley et al. |nj have made it clear that atom-molecule 
coherence can have an important effect on the dynamics of a Bose-Einstein condensate. 
Besides the quantum evaporation process, we thus want to consider this physics in the case 
of the single-pulse experiments as well. We are able to achieve this because very recently 
we have derived an effective quantum field theory that offers a description of the Feshbach 
resonance in terms of an atom-molecule hamiltonian that captures all the relevant two-body 
physics exactly f3~0fl . Apart from a detailed discussion of the condensate collapse, the main 
application of this paper therefore concerns the effect of quantum evaporation in the single- 
pulse experiments and the investigation of the importance of atom-molecule coherence in 
this case. With respect to the latter remark it should be noted that the effect of atom- 
molecule coherence in the case of the condensate collapse will be neglected in the following, 
because the magnetic field is tuned to a far off-resonant value to induce the collapse. 

In view of this the above the paper is organized as follows. In Sec. [TT] we present and discuss 
the theoretical description of quantum evaporation. In Sec. [TTl]we present the applications to 
the condensate collapse, and to the single-pulse experiments with positive scattering length. 
We end in Sec. [TV] with our conclusions. 
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FIG. 1: The scattering length as a function of the magnetic field for ° 5 Rb in the state \f = 2; mj = 
—2). The position of the resonance is indicated by the vertical dashed line. At the horizontal dashed 
line the scattering length vanishes. The dashed-dotted line indicates the background scattering 
length Ob g ~ — 450ao- 



II. QUANTUM EVAPORATION 



In this section we derive the generalized Gross-Pitaevskii equation that includes the 
correction term due to the quantum evaporation of atoms out of the condensate. From this 
result follows a rate equation for the change in the number of atoms in the condensate. 
Sec. H A| is rather technical and may be omitted in a first reading. To facilitate this the final 



rate equation, which is most important for our purposes, is presented in Sec. [I B 



A. Generalized Gross-Pitaevskii equation 



Although the desired rate equation for the number of atoms can also be derived from the 
imaginary-time formalism by means of a Wick rotation to real time, the equation of motion 
for the condensate wave function can not be derived in this manner. Therefore, we use a 
functional formulation of the Schwinger-Keldysh nonequilibrium theory [RS, [36| developed 
in Refs. p7| , 58), ffSfl , from which the equation of motion for the condensate wave function 
follows directly as the equation for the "classical" part of the fluctuating order-parameter 
field. 

Within this formalism, the Wigner probability distribution of the order parameter is 
written as a functional integral over complex fields ?/>*(x, t) and ?/>(x, t). These fields are 
denned on the Keldysh contour C*, which runs from t to t and then back to t^. The 
probability distribution is given by 



P[<f>*,<f>;t] 



iP(x.,t)=4>{x) ( I 

d[rw}ex P \-s[r^] 

V*(x,t)=</>*(x) in 



(2) 



where we absorbed the appropriate initial condition P [</>*, </>; to] m to the measure of the 
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functional integral [39|. The action in the exponent of the integrand is given by 

V>(x',t')>(3) 



S[ifj*,tfj] = dt' dxV*(x',t') 



d h 2 V 2 T 2B (t') 



dt' 2m 



where ^^(x) is the external trapping potential. The interaction is determined by the 
two-body T(ransition) matrix element T 2B (t) = Aira(t)h 2 /m, where a(t) = a(B(t)) is the 
interatomic s-wave scattering length, and m is the mass of one atom. Note that we explicitly 
allowed the scattering length to depend on time. This is experimentally realized by tuning 
the magnetic field near the Feshbach resonance. 

To arrive at an effective action for the condensate wave function, we explicitly separate out 
the condensate part from the field ^/>(x, t). Therefore we write ip^t) = ?/>o(x, t) + ^'(x, t), 
and substitute this into the action in Eq. (Q). In this separation, i/; (x,t) describes the 
condensed part of the gas, whereas ^/>'(x, t) describes the fluctuations. The precise distinc- 
tion between the condensate and the noncondensed part is discussed in detail in Sec fH\ . 
Physically, the idea is that ipo(x.,t) describes the Bose-Einstein condensate and its collective 
modes, whereas ^'(x, t) is associated with the modes not occupied by the condensate. To 
define these two parts consistently, we have to require that they are essentially orthogonal, 
i.e., 

'rfx[^(x,t)V'(x,t) + V / *(x,t)Vo(x,t)] =0. (4) 



This condition ensures that the Jacobian of the transformation of integration variables in 
the functional integral in Eq. (^) is equal to one. In the operator formalism, this condition 
implies that the Bose field operators ?/>'(x, t) and V>'^(x, t) associated with the fluctuations, 
obey the usual Bose commutation relations in the Fock space built upon the states orthogonal 
to Vq( x >*)- 

After this substitution the functional integral becomes 



x J exp {^sy^s, ik>, r , f] + Z j.s'[r , ^ ^ ■ (s) 

Here, we define S'[ipQ, tpo, such that it contains the terms up to quadratic order in 

the fluctuations. We do not retain the terms proportional to (ip') 2 and (ip'*) 2 , since these so- 
called anomalous terms are only needed to describe the collective motion of the condensate. 
Their effect is therefore already included in the action S[ipQ, ip ] for the part of the field that 
describes the condensate. In principle Si nt [il>Q, ipo, V 1 '] contains terms which are either 
of first, third or fourth order in the fluctuations. We neglect, however, the terms of third 
and fourth order, which is known as the Bogoliubov approximation. This approximation is 
valid for the zero-temperature applications under consideration in this paper. Moreover, the 
terms of higher order in the fluctuations describe the interactions among the ejected atoms 
and are expected to be of little importance in determining the ejection rate for the atoms 
expelled from the condensate. 

We write the quadratic action S'[ipQ, ip , ip r *, ip') as 

S'[^l^ ^'\^] = j <W J dx! J df J dx"^*{x',t')hG-\x',t';x",t")^'{x",t"), (6) 
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where we introduced the Green's function G(x, t, x', t') for the fluctuations by means of 

d h 2 v 2 



dt 2m 



-V cxt (x)-2T^(t)|^ (x,t)| 



G(x, t, x', t') = M(x - x')5(t, t'). (7) 



Here, the delta function in the time variables is defined on the Keldysh contour, by means 
of f c t dt'5(t,t') = 1. The part of the action that describes the interactions between the 
condensed and noncondensed parts of the system is, in first instance, given by 



+ V / *(x',t') 



dt' J dx.' 



c* 



8 fi 2 V 2 
<9t' 2m 



^ ext (x')-T 2B (t')|^o(x , ,t / 



V/(x',f) 
V>o(x',t')j- (8) 



It is important to note that this part of the action does not vanish because the condensate 
wave function, as we see in a moment, does not obey the usual Gross- Pit aevskii equation 
once we include the quantum evaporation process. Furthermore, the field ijj'(x.,t) describes 
the high-energy part of the system, and thus has an expansion in terms of the high-energy 
trap states. As a result we are allowed to neglect the terms proportional to the single- 
particle Hamiltonian. The terms with the time derivative vanish because of the orthogonality 
condition in Eq. (|J). The action in Eq. (f|) therefore reduces to 



dt' J dx' [J*(x' } t')^'(x',t') + ^'*(x',t')J(x',t')} 



where we introduced the "current density" 

J(x,t)=T 2B (t)|^ (x,t)| 2 ^o(x,t). 



(9) 



(10) 



The functional integral over the fluctuations in Eq. (^) is a Gaussian integral and can 
be easily performed. This defines the effective action for the condensate on the Keldysh 
contour by means of 



P[<f>*,<f>;t] = J d[r ]d[H^P{j-Sm,H 

x I dW*}d[^] exp V>o, $\ ^) + \S'[^1 Vo, ^ 



d&oWo] exp |-S e f#o^o] j > 
and results in 

— J dlf J dx! J dt" J dx" r(x',t')G(x',t';x",t")J(x",t"). 
Because the Green's function in this effective action is equal to 

iG(x',t';x",t") =Tr{p(t )T et (f(x',t')^(x",t"))} j=Q , 



(11) 



:i2i 



(13) 



where Tc* denotes time ordering along the Keldysh contour and p(to) represents the initial 
density matrix of the gas, this Green's function can be decomposed into its analytic pieces 
by means of 



G(x, t; x', if) = 6(t, t / )G > (x, t; x', if) + 6(t\ t)G<(x, t; x', if), 



(14) 



with 9(t,t') the Heaviside function on the Keldysh contour. The Green's functions G > and 
G < thus correspond to averages of a fixed order of creation and annihilation operators. Note 
that these Green's functions essentially describe the propagation of a noncondensed atom 
in the presence of the mean-field interaction with the condensate. 

To finally derive the generalized Gross-Pitaevskii equation for the condensate wave func- 
tion, we want to separate out the "classical" part of the field "0o( x , t). This is achieved by 
means of the transformation 



V>o(M±) = 0(x,t) ± 



£(M) 



(15) 



Here, </>(x, t) denotes the condensate wave function, whereas the field £(x, t) describes its 
fluctuations. The upper (lower) sign in Eq. ( |I5|) corresponds to the forward (backward) 
branch of the Keldysh contour. When substituting this transformation into the effective 
action in Eq. ([12]), we should in principle only keep terms up to quadratic order in the 
fluctuations, to avoid a double counting of the interactions that we have already taken 
into account. However, to read off the generalized Gross-Pitaevskii equation, including the 
correction terms associated with the quantum evaporation process, it suffices to consider 
only the linear terms in the fluctuations. This is can be understood from the fact that with 
this approximation a functional integration over the fluctuations leads to a constraint for 
0(x, t), which is precisely the "classical" equation of motion that we are interested in. With 
the transformation in Eq. ( JIB] ) we thus project the effective action on the real-time axis 
and read off the equations of motion for 0(x, t) and 0*(x, t) by putting the coefficient of 
the terms linear in £*(x, t) and £(x, t) equal to zero, respectively. After straightforward but 
somewhat tedious algebra this results in 



ih 



gg(M) 

dt 



h 2 v 2 



+ 



2m 

T 2B (t) 

h 



+ V cxt (x) + T 2B (t)|0(x,t)| 2 
dt' f dx'T 2B (t' 



x,t) 

x,t)G (+) (x,t;x',t , )0(x , ,t') 



+<Kx,t)G(-V,t';x,t)0*(x , ,t / ) 



(m) 



(16) 



with the complex conjugate expression for 0*(x, t). We defined the retarded and the ad- 
vanced Green's functions in the usual way by 



G (±) (x, t; x', if) = ±6(±(t - t')) [G>(x, t- x\ if) - G<(x, t- x\ tf)\ . 



(17) 



Note that the Heaviside function in this definition is precisely such that the equation of 
motion for 0(x, t) is causal. 

The time-ordered Feynman diagrams corresponding to the two terms in the general- 
ized Gross-Pitaevskii equation describing the quantum evaporation process are given in 
Figs. |] (a) and (b), respectively. In these diagrams, a condensate atom is denoted by a 
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(a) % (b) 




FIG. 2: These time-ordered Feynmans diagrams correspond to terms associated with quantum 
evaporation in the equation of motion for the condensate wave function 

dashed line. An ingoing dashed line thus corresponds to a factor 0(x, t), and an outgo- 
ing dashed line to a factor <p*(pc,t). The advanced and retarded propagators of the ejected 
atoms are denoted by solid lines. A retarded propagator is denoted by an arrow pointing 
from (x',t') to (x, t), and corresponds to the propagation of a particle. An advanced propa- 
gator is denoted by an reversed arrow, and can be interpreted as the propagation of a "hole" . 
Note that the first term in Eq. flT6|) has an additional factor of two with respect to the second 
term. This is understood from the fact that the corresponding diagram (Fig. ^ (a)) has two 
outgoing lines at (x, t) and therefore contributes twice, whereas the diagram in Fig. |2| (b) 
has only one outgoing line at those coordinates. 



B. Rate equation 

In the previous section we have generalized the Gross-Pitaevskii equation to include the 
quantum evaporation process. With this equation we derive a rate equation for the number 
of atoms in the condensate. It is given by 



dN c (t) d 



dt 

xlm 



r 2B (i)|^(x,i)IV(x,«)G<+»(x, f ;x',i')r 2B («')l^(x',i')IV(x'.i')]. (18) 



The retarded propagator G^ + \x, t; x', t') describes the noncondensed atoms. Note that in 
both this rate equation, and the equation of motion in Eq. ([Tj]) we have taken the limit 
to — > — oo to eliminate initial transient effects. Furthermore, we have made use of the fact 
that G*( + )(x, t; x', t') = G^(x', t'; x, t), which follows from the definition of this Green's 
function. 

The nonmarkovian rate equation in Eq. ([Tj|) describes the evolution of the number of 
condensate atoms due to a change in the scattering length or the condensate wave func- 
tion. Physically the coupling between the condensed and noncondensed parts of the system 
allows the condensate to eject atoms. Since we have neglected the interactions among the 
ejected atoms, this rate equation is applicable only for short times and does not describe the 
rethermalization of the ejected atoms. Because of this approximation the ejection process is 
coherent, and atoms can come back into the condensate on short time scales. To get more 
physical insight in the quantum evaporation process, we discuss the weak-coupling limit, 
where it just corresponds to elastic condensate collisions. 



S 



In the weak-coupling limit, where the mean- field interaction is small compared to the 
energy-level splitting of the external trapping potential, the rate equation for the number of 
atoms in the condensate acquires a familiar form. To see this, we note that in this limit the 
retarded propagator is given by 

G (+) (x, f ; x', = -i8(t - E e-s £n ^Xn(x) X * (x'), (19) 

n^O 

because the single particle states Xn( x ) and energies e n are, in this limit, by definition not 
affected much by the interactions. The sum in this expression is over all the trap states, 
except for the condensate mode, i.e., the one-particle ground state. 
The condensate wave function is given by 

0(x,t)=e-^ i V / iVe(t) Xo (x), (20) 

where fi = eo is the ground state energy, and iV c ^> 1 is the number of atoms in the 
condensate. Using these expressions, the rate equation for the number of atoms in Eq. (|i"8|) 
can be rewritten as 



dN c (t) 
dt 



0-77- A 2 

-^J25(e n -e )\(n\V\0)\ , (21) 

n^O 



which is precisely Fermi's Golden Rule for the rate to scatter out of the initial state |0), 
found from second-order perturbation theory. The sum is over all final states of the form 

(xjxaln) = -j= [xo(xi)Xn(x 2 ) + x n (xi)x (x 2 )] , (22) 

with energy e n + fi. Since we are dealing with identical bosons these states are symmetric. 
This final state thus represents a condensate atom and an ejected atom, whereas the initial 
state with energy 2/i is given by 

(xix 2 |0) = Xo(xi)xo(x 2 ), (23) 

and therefore represents two condensate atoms. The additional factor N% is a result of the 
Bose statistics of the atoms. There is a factor N C (N C — l)/2 « N%/2 for the number of 
condensate atom pairs and an additional factor 1 + N c w N c for the condensate atom that 
is Bose stimulated back into the condensate. Finally, the potential is given by 

U = T 2B 5( Xl -x 2 ). (24) 

We see that, in the weak-coupling limit, the energy-conserving delta function can never 
be obeyed, and thus dN c /dt = 0, as expected in this limit. Another important observation 
is that energy conservation also forbids the ejection of atoms out of a static condensate with 
repulsive interactions, even in the strong-coupling limit. This is because of the fact that 
the energy of a condensate atom, i.e., the chemical potential, is positive in this case, and 
that the energy of the excited states is always larger than the chemical potential. In the 
Thomas- Fermi limit, where the mean-field interaction is much larger than the kinetic energy 
of the atoms, this is illustrated in Fig. [| In this figure, the harmonic trapping potential is 
denoted by the dashed line. For a given chemical potential /i, indicated by the dotted line, 
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FIG. 3: Energy levels of the noncondensate atoms (thin solid lines), in the presence of a condensate 
at a given chemical potential \x (dotted line). The potential felt by the noncondensed atoms is 
indicated by the solid line. 

the noncondensed atoms feel the potential indicated by the solid line. The energy levels 
of the states in this potential are also drawn schematically. Clearly, the energy of the first 
excited state is always larger than the chemical potential, and therefore the ejection of atoms 
out of the condensate by means of elastic collisions is forbidden due to energy conservation. 

This leads to the important conclusion that in order to eject atoms out of the condensate 
by means of the quantum evaporation process, there has to be a strong time dependence 
of either the condensate wave function, or of the scattering length. In the next section, we 
discuss two experimentally relevant examples where this is the case. 



III. APPLICATIONS 



In this section we apply the theory, described in the previous section, to two different 
experimental situations. In the first part, we consider the case where the Feshbach resonance 
is used to investigate the collapse of a condensate 0, |21[|. In the second part we apply our 



theory to the situation where the scattering length is changed very rapidly back and forth 
from zero to a large and positive value [28] . We also discuss the importance of atom- molecule 
coherence in this case. 



A. Condensate collapse 

The use of a Feshbach resonance has made it possible to explore the physics of the 
collapse in two different regimes. First, we can have a collapse in which the dynamics of 
the condensate is mostly determined by the mean-field interactions. This is the case in 
the experiments with 7 Li [|1~J, [16], [H], [Tj|, which has a fixed negative scattering length. In 
these experiments, there is always a large thermal component which feeds the condensate. 
Therefore, we call such a collapse a Type I collapse, in analogy with a Type I supernova, 
which is believed to be the result of an accreting white dwarf that explodes when the 
accumulated mass becomes to large |oT)fl , similar to the 7 Li condensate that grows from 



the thermal cloud and then collapses. The more recent experiments by Donley et al. |2T 
also deal with a Type I collapse, since one starts a collapse from the noninteracting limit, 
where the density is relatively high compared to Bose-Einstein condensates with a repulsive 
self-interaction. 
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Complementary to this regime we can have a second type of collapse where the dynamics 
is mostly determined by the external trapping potential. This regime was considered ex- 
perimentally by Cornish et al. ||. Using the Feshbach resonance of 85 Rb one first makes a 
large, stable and essentially pure condensate, i.e., there is no visible thermal cloud present, 
with repulsive interactions. Then, one suddenly switches the interactions from repulsive to 
slightly attractive, and watches the subsequent collapse. Because the collapse in this case 
starts at low density due to the initially large and repulsive interatomic interactions and the 
magnitude of the final negative scattering length is much smaller than the initial positive 
one, its dynamics is mostly determined by the trapping potential. We call such a collapse 
a Type II collapse, in analogy with a Type II supernova. Such a supernova is the fate of a 
massive star analogous to the large condensate with repulsive interactions. 

In the experiments on the Type II collapse it was found that after such a collapse the 
number of atoms in the condensate is of the same order as the maximum number of atoms al- 
lowed to have a metastable condensate j|(J . This suggests that during such a collapse enough 
atoms are ejected from the condensate, so that the condensate becomes metastable. The 
most important theoretical task is therefore to identify the physical mechanism responsible 
for this ejection. We have recently argued that the quantum evaporation process provides 



an explanation for these experiments p9fl . Moreover, we will explicitly show here that the 
decay solely by means of three-body recombination does not explain the experimental results 
for the Type II collapse. 

Interestingly, the metastable condensate that is the result of a Type II collapse is pre- 
cisely the starting point for a Type I collapse, making this last collapse a much more violent 
phenomenon. We restrict ourselves here to the description of the Type II collapse, because 
a proper treatment of the Type I collapse requires a full numerical solution of the general- 
ized Gross-Pitaevskii equation that includes the nonlocal correction term due to quantum 
evaporation. This is beyond the scope of the present paper, and we use a simpler variational 
approach here that we believe is appropriate for a Type II-collapse event. 



1. Gaussian approximation and semiclassical retarded propagator 

In principle, we must now numerically solve the generalized Gross-Pitaevskii equation in 
Eq. that includes nonlocal terms. However, in order to gain physical insight, we will 
make several approximations to reduce the solution of this equation to a numerically more 
tractable problem. First of all, we use a gaussian variational approach to the condensate 
wave function pT| . More precisely, we assume the condensate wave function to be of the 
form 




') = ^X(/K'^"n(^77T) xox P f--- T ^(:i-,^ Z ^pj| . (25) 

In the limit of a small number of condensate atoms N c (t) this ansatz becomes an exact 
solution of the Gross-Pitaevskii equation for a harmonic external trapping potential, and 
therefore a good description of the condensate after the collapse. It is, however, also known 
that a gaussian ansatz gives good results on the frequencies of the collective modes, even 



in the Thomas- Fermi regime ||41|| . Moreover, since we consider only the Type II collapse 
and therefore by definition assume that the trapping potential is more important than the 
mean-field interactions, we expect the gaussian ansatz to give physically reasonable results 
at all times. 
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In the case of the ordinary Gross-Pitaevskii equation the variational parameters obey 
Newton's equation of motion [12 



m^M = Sv( q (t);N c (t)), (26) 



dt 2 dq 



j 



with a potential given by 

where the frequencies uij are the frequencies of the harmonic external trapping potential 
in the three spatial directions. We have recently extended the above variational calculus 
also to the case of a nonlinear Schrodinger equation including an imaginary term due to 
the presence of a thermal cloud, and found that the equation of motion for the variational 



parameters in principle contains a damping term ||18| . However, we are interested in the 
zero-temperature situation, in which this damping term is negligible. This means that for 
a description of the collapse, we only have to couple Eq. ( p6|) to the rate equation for the 
number of atoms in Eq. flT8|). 

To determine this rate equation in detail, we need an expression for the propagator of 
the ejected atoms. Anticipating that the energy of the ejected atoms will be much higher 
than the energy of a condensate atom, we take for the one-particle Green's function of the 
ejected atoms the semiclassical approximation, 

G«(x,t;x',t') = -i9(t-t') f ^e-HM^Xt-tVMx-x^ (2g) 

where we introduced the center of mass coordinate R = (x + x')/2, and used T = (t + f)/2 
for notational convenience. The energy of the ejected atoms is, in this approximation, given 
by 

e(k, R, T) = + V ext (R) + 2T 2B (T)|0(R, T)| 2 + e x . (29) 
2m 

This energy is measured from ei, which is the first excited level in the trap. This energy 
thus determines the cut-off between the condensate and the ejected atoms. Although this 
cut-off is important to make the distinction between condensate and noncondensate atoms, 
we find that our numerical results are not very sensitive to its precise value. 

Since the most important contribution to the rate equation comes from the center of the 
trap, we expect that we can, to a good approximation, take for the energy of an ejected 
atom the value e(k, 0,T) = e(k, T). We use this value for the energy of an ejected atom 
from now on. This approximation is very convenient, since we can now perform the spatial 
integrals in the rate equation, as well as the integral over the momenta in the expression for 
the retarded propagator analytically. The final result is given in Appendix [A[ 

In the ansatz in Eq. (^) we have included a global phase #o(t). This is important, since 
we have already seen in the previous section that this global phase determines the energy 
of a condensate atom in equilibrium. We can determine this global phase by insertion of 
the ansatz in the Gross-Pitaevskii equation including the correction term in Eq. fli~6|) , and 
separating out the real part. Neglecting the contribution of the quantum evaporation process 
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to the mean-field energy, it is up to an irrelevant constant given by ]T8 



9 (t) 




dr 



dV(q(T),N c (r)) 
dN r 



El -2 
-A m( lj 



Note that in equilibrium this results in a phase factor e~ %fJ ' t ' n for the condensate wave func- 
tion, with fi = dV(q, N c )/dN c . This is clearly the correct expression for the chemical 
potential. 

With the above approximations we have reduced the difficult problem of solving a partial 
differential equation with a nonlocal correction term to the problem of solving the ordinary 
differential equations for the variational parameters in Eq. fl26|) coupled to the rate equation 
for the change in the number of atoms in Eq. (18) . This rate equation is now only nonlocal 



in time, since all the integrals can be done analytically for the gaussian wave function and 
the semiclassical retarded propagator we are considering here. In the next section we will 
present our results obtained by numerically solving these equations. 



2. Results 

We have performed numerical simulations for the experimental conditions of Cornish 
et al. [|9], |4(| with 85 Rb. The frequencies of the external trapping potential are equal 
to ijj r /2 , K = 17 .4 Hz and u z /2tt = 6.8 Hz, in the radial and axial direction, respectively. 
One starts with a condensate consisting of N c (0) ~ 4000 atoms and a scattering length of 
a(-Bj) = 2500ao, where B x is the magnetic field at the initial time. Initially, there is no visible 
thermal cloud present so the system is approximately at zero temperature. One then ramps 
the magnetic field linearly in a time At from its initial value B\ to the final value Bf, chosen 
such that a(Bf) = — 60ao- The ramp time is taken equal to At = 0.5 ms. 

To investigate the importance of three-body recombination during the collapse, we have 
included it in our simulations in addition to the quantum evaporation process. This amounts 
to including a term 

-|ffl0(x,t)|V(x,t) 

on the right hand side of Eq. (|16"D, from which the contribution to the rate equation for the 
number of atoms can be easily found. Although there are several predictions for the normal- 
component rate constant K3 |f42| , fE| [EJ, its behavior as a function of the magnetic 



field is unknown near the Feshbach resonance and precise experimental data is unavailable 
Following Saito and Ueda | ZEfl , we take K 3 = 1.2 x 10~ 27 cm 6 /s. With this value 



these authors have been able to explain some of the results of the experiment on the Type 
I collapse by Donley et al. |2l|]. Since the final value of the magnetic field is in the same 



range for both experiments, we expect the three-body recombination rate constant to be of 
the same order of magnitude. 

In Figs. f| and |5| we present the results of our simulations. The widths of the condensate 
in the radial and axial direction, i.e., the variational parameters q r and q z , are shown as a 
function of time in Fig. |] (a) and (b), respectively. The number of atoms as a function of 
time is displayed in Fig. |5|. The origin of the time axis is chosen such that the scattering 
length is equal to zero at t — 0. The solid line corresponds to a simulation that includes both 
the quantum evaporation process and three-body recombination. The condensate collapses 
first in the radial direction in a time n/(2u) r ) ~ 14 ms and during the last part of this 
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FIG. 4: (a) Radial and (b) axial width of the condensate as a function of time during a Type II 
collapse. At the origin of the time axis the scattering length vanishes as it is changed from large and 
postive to negative. The solid line corresponds to a simulation where both quantum evaporation 
and three-body recombination are included. The dashed line corresponds to a simulation including 
only three-body recombination. The experimental points are also shown and taken from Ref. [40 1. 



radial collapse the condensate ejects a large fraction of its atoms, by means of the quantum 
evaporation process. In our simulations, we find that the three-body recombination hardly 
contributes. After a time ir/(2u z ) ~ 36 ms the axial width of the condensate reaches its 
minimum, resulting again in a slight increase of the ejection rate. Note that these time 
scales are expected, since the dynamics is determined by the external trapping potential. 
More remarkably, for a very short time the number of atoms increases with time during 
the collapse. In principle, this can occur, because the quantum evaporation process is 
coherent, and can thus describe multimode Rabi oscillations between the condensate and 
the noncondensed part of the system, as we shall see in much more detail in the next section. 

The simulation that includes the quantum evaporation process shows quantitative agree- 
ment with the experimental results [40|, which are also shown in the Figs. f| and |5|. The 
disagreement for small values of the widths of the condensate is a result of the fact that the 
experimental resolution for the condensate size is about 4 jum [40]. 

We have also performed a simulation that includes only three-body recombination and 
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FIG. 5: Number of atoms in the condensate as a function of time during a Type II collapse. The 
solid line corresponds to a simulation where both quantum evaporation and three-body recom- 
bination are included. The dashed line corresponds to a simulation including only three-body 
recombination. The experimental points are also shown and taken from Ref. pOp. 



no quantum evaporation. The results for this simulation are shown in Figs. f| and [| by 
the dashed lines. The fact that the minima in the widths of the condensate in Fig. |] are 
lower than the results including quantum evaporation, indicates that the condensate density 
has to be relatively high for recombination of atoms to occur. Once such a high density is 
reached, the ejection of atoms occurs very fast, resulting in a staircase- like pattern for the 
number of atoms as a function of time, which is clearly not visible in the experimental data. 
At t ~ 38 ms the condensate decays so fast, that the radial direction becomes stable again, 
resulting in a large increase of the radial width. 



Finally we make some remarks about the properties of the ejected atoms. In Ref. |29 
we have calculated the distribution of the kinetic energy of the ejected atoms, as well as 
the angular distribution. We have also performed preliminary calculations of the average 
kinetic energy emitted in the radial and axial directions, and have found that the energy 
distribution of the ejected atoms is anisotropic. Donley et al. [21] have measured the 
angular and radial temperatures of the emitted atoms and have indeed found an anisotropic 
distribution of the energies. However, our calculations are done for the Type II collapse 
and our approximations are especially suited for this case, whereas these experiments deal 
with a Type I collapse. Therefore, we do not directly compare our results with the available 
experimental data of Donley et al. [BTj]. Moreover, to improve the experimental resolution 
limit, one expands the gas at the end of each destructive measurement by an increase in the 
scattering length. Therefore, the resulting mean-field energy and rethermalization effects 
may play an important role in determining the energy of the ejected atoms, which is not 
included in our calculations of the energy distribution function in Ref. |29 



In conclusion, we have shown in this subsection that the simulations of the Type II col- 
lapse that include quantum evaporation show quantitative agreement with the experimental 
results. The most important feature of these experimental results is the fact that the con- 
densate starts to eject atoms almost immediately after the initiation of the collapse. We 
have also shown that solely three-body recombination does not account for this rapid onset 
of the loss of atoms. At this point it is important to notice that this conclusion also holds 
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FIG. 6: The scattering length as a function of time in the experiments by Claussen et al. |28|| , 
and the corresponding magnetic field (inset). One increases the scattering length in a time t T i se by 
means of a linear ramp in the magnetic field, and holds the magnetic field for a time thoid 5 before 
ramping back. The position of the Feshbach resonance is indicated by the dashed line. 



if we numerically solve the generalized Gross-Pitaevskii equation, without approximations. 
The reason for this is that the gaussian ansatz used here is certainly appropriate for the 
first part of the collapse, when the dynamics is not yet very violent. This is borne out by 
numerical simulations of the Gross-Pitaevskii equation, which have also shown that when 
the highest densities are reached during a Type I collapse, high-density "spikes" can form 



on the profile of the wave function [23, G4L Ed, E6l K7j. These "spikes" are not included 



in our ansatz for the condensate wave function and therefore our approximations might be 
less appropriate for the highest densities reached during the collapse. However, our results 
suggest that such high-density "spikes" may well never occur if one includes the effects of 
the quantum evaporation process. 



B. Multimode Rabi oscillations 

Apart from the negative scattering length regime, the experimental control over the 
interatomic interactions has also made it possible to explore the regime where the interaction 
is large and positive. To this end, Claussen et al. [E3[] have conducted an experiment where 
the magnetic field undergoes a trapezoidal pulse in time, resulting in a quick jump in the 
scattering length towards a large positive value. In detail, one ramps the magnetic field 
linearly in a time t r i se from its value in the noninteracting limit (B « 166.5 G), to a value 
where the scattering is of the order of a few thousand Bohr radii. The magnetic field is kept 
at this value for a time thoid before ramping back to the initial value again within the same 
t rise . The scattering length as a function of time for a typical pulse is shown in Fig. |6] and 
the inset shows the corresponding magnetic field. The rise time t Tise and the hold time ihoid 
are typically of the order of 10 — 100 /is. 

In this experiment, one observes particle loss from the condensate as a function of the both 
the rise time and the hold time, accompanied by a "burst" of atoms from the condensate. 
The temperature of the burst atoms is of the same order as in the case of the experiments 
on the collapse, i.e., about 150 nK. Most importantly, the amount of atoms lost from the 
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condensate decreases with increasing t r ise over some interval, which can never be the case 
for atom loss characterized by a rate constant, such as three-body recombination. Also note 
that the time scales of the pulse are very small compared to the time scales set by the inverse 
frequencies of the external trapping potential. This means that the density profile of the 
condensate hardly changes its shape during the pulse, since the collective modes which alter 
the density profile have frequencies on the same order of magnitude. 



1. Retarded propagator 

The single-pulse experiments start in the noninteracting limit where the density profile 
of the condensate has the shape of a gaussian, and the pulse in the magnetic field is so fast 
that the condensate wave function hardly changes due to the interactions. This makes it 
convenient to expand the propagator for the ejected atoms in the excited harmonic oscillator 
eigenstates, since the condensate part is then left out most easily. The fact that the con- 
densate density profile almost remains the shape of the ground state of the trap makes this 
problem easier to deal with theoretically and makes it therefore worthwhile to determine the 
propagator for the ejected atoms as accurate as possible. With the Type II-collapse problem 
of the previous subsection this objective is much more difficult since then the condensate 
wave function changes its shape considerably during the collapse. This means that at each 
time a different number of trap states has to be included in the wave function, making 
the cut-off between condensate and noncondensed atoms strongly time- dependent. There- 
fore we have applied several approximations in that case. Even though the density profile 
of the condensate does not change much, the phase of the condensate wave function does 
change very fast in the single-pulse case. Therefore we use for the phase of condensate wave 
function the phase of the gaussian ansatz given in Eq. (p5|). Since this wave function also 
describes the low-lying collective modes of the condensate, we have to exclude these from 
the propagator of the noncondensed atoms. 

Denoting the eigenstates of the single-particle hamiltonian again by Xn( x ) and the cor- 
responding eigenvalues by e n , we have for the propagator of the ejected atoms 



GWfx,^/) = -i9(t-t'] 



n,n' 



«(t,t / )Xn(x)x^(x')e- 



4(e n t-e n ,t') 



(31) 



The prime denotes summation over all the excited trap states not contained in the condensate 
wave function. The equation of motion for the expansion coefficients a njn '(t, t') can be found 
by inserting Eq. (|3l| ) into the equation of motion for the Green's function given in Eq. (^). 
However, it is easier to realize that the Green's function is in Eq. ( |I3"D shown to be related to 
the expectation value of the product of two Heisenberg annihilation and creation operators 
for the noncondensed atoms. The equation of motion for the annihilation operator of interest 
is given by 



d h 2 v 2 



2m 



-y cxt (x)-2T^(t M x,t)| 



f(x,t) = 0, 



(32) 



with the hermitian conjugate expression for the creation operator. We solve this equation 
by expanding the annihilation operator as 



^'(X,t) = X!^m(x,t)V4 



(33) 
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FIG. 7: The fraction of atoms in the condensate as a function of time for a calculation that 
includes only quantum evaporation. The initial number of condensate atoms is iV c (0) = 16500. 
The rise time i r i se = 12.5 fis and the hold time is ihoid = 200 [as. The scattering length is equal to 
a = 2000ao during hold. 



where the Schrodinger operator ip' m annihilates an atom in the harmonic-oscillator state with 
quantum number m. We then expand also the functions m (x, t) in trap states by means of 

m (x ) t)=5: / C(t)e-^"'xn(x), (34) 

n 

and determine the equation of motion for the coefficients c™(t) from Eq. (p2|). This results 
in 

d ^m_ = -|T 2B (t)EV n , n ,(t)C(t)e-i^-) t) (35) 

with matrix elements given by 

V^it) = J rfx X ;(x)|0(x,t)| 2 Xn ,(x), (36) 



which depend on time through the variational parameters in the gaussian ansatz in Eq. ( |25|) 
and the number of condensate atoms. These matrix elements can be calculated analytically 
and the result is given in Appendix. [B]. The advantage of the above approach is that we do 
not have to solve the equation for the derivative of a ntT1 '(t,t') with respect to t' separately. 

Putting the results together, we find for the coefficients in the expansion of the Green's 
function the expression 

an,„'(M') = E'CWC(0- (37) 

m 

With this Green's function we have performed simulations of the single-pulse experiments 
by Claussen et al. [p^j, of which the results are presented in the next section. 



2. Results 

We perform our calculations for the parameters of the experiment by Claussen et al. 



28] . In particular, the frequencies of the external trapping potential are the same as in the 
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FIG. 8: Fraction of atoms in the condensate as a function of the hold time. The rise time is kept 
fixed at the value t r i se = 12.5 /is. The initial number of condensate atoms is taken -/V c (0) = 6100 
(solid line) and N c (0) = 16500 (dotted line). The scattering length is equal to a = 2000ao during 
hold, (a) The result of the calculation that includes only quantum evaporation, (b) The calculation 
that includes both quantum evaporation and three-body recombination. The experimental points 
are taken from Ref. pg|]. 



previous section. Fig. [7] shows the fraction of atoms in the condensate as a function of time, 
for a pulse such that t r i se = 12.5 /xs and thoid = 200 /xs. The magnetic field during the hold is 
B = 156.9 G, which corresponds to a scattering length of a = 2000ao- The simulation shows 
that once the scattering length nearly takes on its largest value and the coupling between 
the condensate and the excited states is therefore largest, the condensate starts ejecting 
atoms. Part of these atoms then oscillate back and forth between the condensate and the 
excited states. The curve in Fig. [7| clearly contains several frequencies since we are dealing 
with several excited states and thus a multimode Rabi oscillation. At the end of the pulse 
the rate, i.e., the slope of the curve, becomes equal to zero because the coupling between 
the condensate and the excited states becomes equal to zero at the end of the pulse, where 
the scattering length is equal to zero. 

To compare our results with the available experimental data we calculate the number of 
atoms as a function of the hold time thoid and the rise time t rise . Fig. (a) shows the result 
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FIG. 9: Fraction of atoms in the condensate as a function of the rise time for a calculation that 
includes only quantum evaporation. The hold time is kept fixed at thoid = 1 jus (solid line), ihoid = 5 
/jls (dashed line), and ihoid = 15 fjs (dotted line). The scattering length is equal to a = 2000ao 
during hold. 



of a calculation of the fraction of condensate atoms as a function of thoid, with t r i se = 12.5 
/is fixed. The calculation is done for two different initial numbers of condensate atoms. The 
solid line displays the result for N c (0) = 6100 and the dashed line for iV c (0) = 16500. Notice 
that the latter initially has a larger slope because the effective Rabi coupling between the 
condensate and the excited states is larger in this case. This is because of the fact that it is 
proportional to the condensate density. 

The results of the simulation that includes only quantum evaporation, shown in Fig. |8| (a), 
show an oscillation in the fraction of condensate atoms as a function of the hold time. This 
oscillation is not observed in experiment, because of the fact that three-body recombination 
plays an important role in this case since it becomes large near the resonance [46|. Therefore, 
we also want to perform a calculation that includes both quantum evaporation and three- 
body recombination. However, the magnetic-field dependence of the rate constant for this 
process is unknown. Nevertheless, we are able to make progress by realizing that the hold 
time is generally larger than the rise time for the experimental points shown in Fig. |8|. Since 
experimentally three-body recombination is known to increase by orders of magnitude near 
the resonance |fi6fl , the contribution of three-body recombination will be most important 



during hold, where the magnetic field is closest to the resonance. This suggests that we only 
need to include it during hold. Note that for this approximation to be valid it is essential that 
the rise time is shorter that the hold time. If the rise time is larger than the hold time the 
magnetic field dependence of the three-body recombination rate constant is of importance, 
since the magnetic field is then time dependent for almost the entire pulse. 

Fig. |S] (b) shows the result for a calculation that includes both quantum evaporation 
and three-body recombination, with a rate constant K 3 = 3 x 10~ 23 cm 6 /s during hold. 
This value for the normal- component rate constant agrees with the order of magnitude 



of the experimental data f46fl . This calculation shows good quantitative agreement with 
experiment for both initial numbers of condensate atoms. 

Finally, we have calculated the number of atoms in the condensate as a function of the 
rise time. The result of this calculation is shown in Fig. || for various hold times. The solid 
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line corresponds to thoid = 1 £*s. The dashed and dotted line correspond to hold times of 
5 /is and 15 /is, respectively. The scattering length during hold is equal to a = 2000ao for 
this simulation. For all the results in Fig. ^| the rise time of the pulse is larger than the hold 
time. This means that the magnetic-field dependence of the three-body recombination rate 
constant is very important in this case, since the magnetic field is varying most of the time. 
Fitting the dependence to the experiments is difficult due to the long times taken by the 
numerical computations. Therefore we refrain from including three-body recombination in 
these simulations. Nevertheless, there is agreement with the experimental results regarding 
several aspects of our results. First, we find that the number of atoms increases with the 
rise time over some interval. This was also found in the experiment by Claussen et al. p8 



Note that this fact can not be explained by any loss process characterized by a rate constant 
because the amount of atom loss will then always be larger with longer times. Second, the 
minima of the curves in Fig. ^ shift to lower values of t r i se with longer hold times. This 
was also observed in the experiment by Claussen et al. [28]. These minima also occur on 
approximately the experimental values of t r - lse . The fact that in the experiment the minima 
become lower with increasing hold time can be explained by three-body recombination. 

In conclusion, we have applied the generalized Gross-Pitaevskii equation to the recent 
single-pulse experiments by Claussen et al. 



281. We have shown that the number of atoms 



increases with time over some ranges. This can not be understood in terms of conven- 
tional loss processes such as three-body recombination or dipolar decay. However, to obtain 
agreement with the available experimental data we had to include three-body recombination 
in our simulations. Due to the fact that the magnetic field dependence of this process is 
completely unknown we are not able to make a fit to experiment in all situations. 

C. Atom- molecule coherence 



Recent experimental and theoretical work has shown that atom-molecule coherence is of 

Therefore 



importance in the case of a double pulse in the magnetic field pi], [32], |34 



we 



may expect it to have an important effect in the case of the single-pulse experiments as well. 
To make the discussion of these experiments more complete, we investigate the role of the 
molecules by means of a quantum field theory that we derived recently pOfl . This theory 
incorporates the correct molecular binding energy and scattering properties of the atoms 
at the quantum level by using coupling constants that are dressed by ladder diagrams and 
by including the molecular self-energy. Introducing Heisenberg operators ip a and ip m that 
annihilate an atom and a bare molecule, respectively, the hamiltonian for the atom-molecule 
system reads 



ih 



ih 



dip m (x,t) 
dt 



dt 



2V72 



rv 

4m 

ft 2 V 2 
2m 



+ 6(B(t))-g 



2 m 3 / 2 Id h 2 V 2 
2nh 3 V dt Am 



+ T b 2B ^(x,t)^ a (x,t) 



^ m (x,t) +gipl(x,t) , 



^ a (x,t) + 2#t(x,t)</> m (x,t) 



(38) 



Here, g = h^Tra^gABAfi/m is the atom-molecule coupling constant and 
6(B) = Afi(B(t) — Bo) denotes the detuning, i.e., the energy difference between two atoms 
and the bare molecule. It is determined by the difference in magnetic moment between the 
atoms and the bare molecule, which in the case of 85 Rb is equal to A/x ~ — 2.2/xb |52 . 
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FIG. 10: Fraction of atoms converted in the condensate as a function of the hold time. The rise 
time is kept fixed at the value i r i se = 12.5 [is. The initial number of condensate atoms is taken 
N c (0) = 6100 (solid line) and N c (0) = 16500 (dashed line). The scattering length is equal to 
a = 2000ao during hold, (a) The result of the calculation that includes only the coupling of the 
atomic condensate to the molecular field, (b) The calculation that includes both atom-molecule 
coupling and three-body recombination. The experimental points are taken from Ref. [28|. 



At first glance the term proportional to yihd/dt + fr 2 V 2 /(4m) may appear unexpected. 
It corresponds to the imaginary part of the self energy of the bare molecule which arises 
physically from the fact that the molecular state interacts with the two-atom continuum. 
This affects both the wave function of the dressed molecule and its binding energy. By 
determining the pole of the molecular propagator for negative detuning, the latter can be 
shown to be given by |J0| 



e m (B) = 5(B) 



4 3 

g m 
8rf 



\ 



167T 2 /i 6 

i - ^-s(b) 



(39) 



which reduces to e m (B) = —h 2 /[m(a(B)) 2 ] for values of the magnetic field close to the 
resonance. Due to the coupling with the continuum of atoms, i.e., the open channel of the 
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Feshbach problem, the molecular state is strongly affected and is given by 



i f c?k 

\Xm, dressed) = ^JZ(B)\x m ; bare) + / jz -^C(k)|k, -k; open) 



(40) 



(2tt)* 

where the coefficients C(k) are normalized as / <ik|C(k)| 2 /(27r) 3 = 1 — Z{B). It contains 

with an amplitude Jz(B) the bare molecular state |x m ;bare). Moreover, because of the 
coupling to the two-atom continuum the molecule acquires a nonzero component in the open 
channel. The wave function renormalization factor Z(B) is given by [3(| 



Z(B) 



1 



g*m 3 / 2 /(iirhy\e m (B)\) 



(41) 



which approaches one for values of the magnetic field far off resonance, where the dressed 
molecular state reduces to the bare molecular state, as expected. However, for values of 
the magnetic field close to the resonance, it is much smaller than one. In particular, for 
the case of the single-pulse experiments we always have that Z(B) <C 1, which implies that 
the magnetic moment of the dressed molecule is in very good approximation equal to twice 
the magnetic moment of an atom. For magnetically trapped atoms this implies that the 
dressed molecule is subject to twice the trapping potential for the atoms. With respect to 
this remark it is important to note that the result of the calculations of Kokkelmans and 
Holland p2] for the density of the molecular condensate should be multiplied by a factor 
1/Z(B) ^> 1 to obtain the density of real dressed molecules, since these authors calculate 
the expectation value of the bare molecular field operator (ij; m (x, t)). 

To bring out the physics of Eq. (|38|) more clearly, we introduce the operator ij)' m = 
4> m / J Z(B) that creates a dressed molecule, i.e., a molecule with an internal state as in 
Eq. fl40D . Since we intend to consider the situation where initially all atoms are in the 



atomic condensate, we are allowed to make a mean-field approximation for the atomic field 
operator and consider only its expectation value. There are however no molecules present 
at the initial time and this requires a quantum treatment of the molecular field operators. 
The resulting equations for the atomic condensate wave function coupled to the dressed 
molecular field reads for the experimental conditions of interest 



ill 



(X,t) 



ih 



dt 

3&(x,t) 

dt 



h 2 V 2 

2m 

h 2 v 2 

4m 



+ V- xt (x)+T b f|0 a (x,t)| 



a (x,t) + 2WZ(t)0:(x,tK4(x,t) 



+ 2^ ext (x) + e m (t) 



<(x,t) + W£(WCM) , 



(42) 



where a = (ip a ), Z(t) = Z(B(t)), and e m (t) = e m (B(t)). In the derivation of the above 
coupled equations we have assumed that we are allowed to make an adiabatic approximation 
for the renormalization factor Z{B) and that we can evaluate it at every time at the magnetic 
field B(t). In principle there are retardation effects due to the fact that the dressed molecular 
state does not change instantaneously. It turns out that these effects can be neglected if 



h 



d\nZ{B{t)) 
dt 



< K(B(t))\ 



(43) 



23 




100 

rise time (us) 



200 



FIG. 11: Fraction of atoms in the atomic condensate for the situation where the atomic condensate 
is coupled to the molecular field. Initially there are iV c (0) = 16500 atoms and no molecules. The 
magnetic field is such that a = 2000ao during the hold, (a) Fraction of atoms as a function of the 
real time for t r j se = 200 /is (solid line) and t r \ se = 100 fis (dashed line). The hold time is equal to 
^hoid = 1 M s f° r both pulses, (b) Fraction of atoms as a function of the rise time for different hold 
times of ihoid = 1 /xs (solid line), ihoid = 5 fis (dashed line) and ihoid = 15 A*s (dotted line). 



which is fulfilled for almost the entire duration of most of the pulses in the experiments of 
Claussen et al. We come back to this point in the discussion at the end of the paper. In 
principle, the coupling between the two-atom continuum and the molecule also contains an 
incoherent part corresponding to the rogue dissociation process considered by Mackie et al. 
33fl . The rate for this process will be small, however, under the condition given in Eq. 



Moreover, the mean-field effects of the condensate on the thermal atoms will suppress this 
process even further. It can, in principle, be included straightforwardly and will take the 
form of a dissipation term in the equation for the molecular operator. 

We solve the equations for the atomic condensate wave function coupled to the dressed 
molecular field by using for the condensate wave function again the gaussian ansatz in 
Eq. (|25|), and by expanding the dressed molecular annihilation operator in harmonic oscil- 
lator eigenstates, similar to the expansion in Eq. (^). As initial condition we assume that 
at t — only condensed atoms are present. The results of our calculations are shown in 
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Figs. ITO and |TT 



The calculations presented in Fig 10 are performed for the same experimental conditions 
as in Fig. || This result clearly shows that a large fraction of atoms is coherently converted 
into molecules as a result of the fast ramp in the magnetic field and that these oscillate 
back and forth between the atomic condensate and the molecular states. Due to the fact 
that the conversion of the atoms to the molecular states is coherent, the operator ^(x, t) 
acquires a nonzero expectation value. Fig |TU| (b) shows the results of simulations where also 
three-body decay is taken into account in the same manner as in the previous section. The 
normal component rate constant is taken equal to K3 = 3 x 10~ 23 cm 6 /s. Interestingly, 
the initial decay without three-body recombination is already larger than the experimental 
data and by adding three-body recombination it is therefore impossible to make a fit to 
the experimental data. This is possibly the result of neglecting the retardation effects of 
the renormalization factor Z(B(t)) and the rogue-dissociation process, since the condition 
in Eq. (|43"D is violated for a significant fraction of the total duration of the pulse in this 



case. For the simulations presented in Fig. |TT] this condition is violated only for a very small 
fraction of the total duration of the pulse for rise times larger that t rise ~ 50 /is and is not 
violated at all for t rise > 150 /is. Note that the effect of retardation and rogue dissociation 
lead to decoherence, which means that our calculations give an upper bound on the amount 
of molecules that are coherent with the atoms. 

Fig. [ll| (a) shows the result of two calculations for different rise time as a function of the 
total time. As expected, the number of atoms in the atomic condensate first oscillates with 
large frequency since the dressed molecular binding energy is large here. As the magnetic 
field approaches values closer to the resonance, the frequency decreases. From Fig. [11] (a) it 
is clear that only the largest frequency, which also has the largest amplitude since the gas is 
then closest to resonance, gives a significant contribution to the frequency observed in the 
number of condensate atoms as a function of the rise time, because the larger frequencies 
with smaller amplitude average out. However, these oscillations are not observed in the 
experimental data of Ref. Introducing three-body recombination to fit the theory to 

experiment is impossible with pulses having relatively long rise times, for the same reasons 
as in the previous section. Nevertheless, the amplitude of the oscillations are in this case, 
except for the longest hold time, comparable to that of the simulations where only quantum 
evaporation is included. This implies that for a thorough treatment of the single-pulse 
experiments both atom-molecule coherence and quantum evaporation should be included. 
This is beyond the scope of the present paper but work in this direction is in progress. 



IV. CONCLUSIONS 

We have put forward a generalized Gross-Pitaesvkii equation that includes nonlocal terms 
which describe the quantum evaporation of the Bose-Einstein condensate. We have applied 
this equation to two experimental situations which make use of a Feshbach resonance to 
alter the interaction properties of the atoms. First we have considered the case where 
the condensate undergoes a Type II collapse whose dynamics is mainly determined by the 
external trapping potential and have found good quantitative agreement with experiment. 
Second we have considered the recent single-pulse experiments p8 |. In general we have 



also found agreement with experiment in this case, keeping in mind that the magnetic-field 
dependence of the three-body recombination rate constant is completely unknown. The 
latter is the first serious complication in the theoretical analysis. Apart from considering 
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quantum evaporation we have also considered the role of atom-molecule coherence in the 
single-pulse experiments, by means of an adiabatic approximation to our effective quantum 
field theory for the description of Feshbach resonances JJIJ . In first instance atom- molecule 
coherence appears to be an important effect. However, the second theoretical complication 
is that the adiabatic approximation in general overestimates the effect and does not take 
into account rogue dissociation Including this process damps out the Rabi oscillations 
between the atoms and molecules and leads to the production of energetic atoms that may 
contribute to the experimentally observed bursts |TT]. Due to these two complications, a 
completely satisfying quantitative description of these experiments is still lacking. It should 
be mentioned that our calculations take into account the inhomogeneity of the trapped gas 
exactly and not in local-density approximation. In addition, we do not make a single-mode 
approximation either for the atomic condensate or the dressed molecules. In future work we 
intend to consider quantum evaporation, rogue dissociation, and three-body recombination 
simultaneously to obtain more insight into these intriguing JILA experiments. 
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APPENDIX A: RATE EQUATION FOR TYPE II COLLAPSE 

With the gaussian ansatz for the condensate wave function and the semiclassical propa- 
gator for the ejected atoms, the final rate equation for the change in the number of atoms 
is given by 



dN c (t) 
dt 



Ui^a{t)Nl' 2 {t) / dt' a(t')N^ 2 (t')e-^- e "^-T \<Ko,(t+t>)/2)W-t) 



{nmq r (t)q r (t') [iq r (t)q r (t') ((t - t')q r (t) - q r (t)) q r (t') - q r (t')q r (t)] m 2 
+3mh [q 2 \t) + (t' - t){q r {t)q r {t) - q r (t')q r (t')) + q 2 (t')] + 9ih 2 (t' - t)} 



LU, 



\ 



q z {t)q z {t') 3 



imq z (t')q z (t')' 
k 



x 



-K: 



mq z (t)q 3 z (t) 
^ 3ihmq z (t)q z (t) 



+ q 2 z (t) + 3ih 



t - v 



3ih + mq z (t')q z (t') m 

t3 3^37^(^(0' 

Here, q r (t) and q z (t) denote the radial and axial width of the condensate, respectively. 



(Al) 
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APPENDIX B: MATRIX ELEMENTS 



In this Appendix we calculate the matrix elements V njn i(t) in Eq. (|36|). Because of the 
fact that we are dealing with a trapping potential that is symmetric around the z-axis, 
the excited states factorize into a radial and an axial part. It is convenient to characterize 
the radial part of the excited states by the quantum numbers (n r , m), where n r counts the 
number of radial nodes in the wave function and m is the quantum number corresponding 
to the projection of the angular momentum on the z-axis. The third quantum number n z 
counts the number of nodes in the axial direction. In cylindrical coordinates these states 
are given by 



Xn M (r,0,z) oc e-^^lr^Fti-nrAml + l^r/lr^e 



imd 



xH nz (z/l z )e 



-*7(2«) 



(Bl) 



where U = yh/(muji). The Hermite polynomials are denoted by H n (x) and the confluent 
hypergeometric function is denoted by iFi(p,q,x). The overlap integral with two functions 
of the form as in Eq. ( |B"T| ) with a gaussian of arbitrary width is, to the best of our knowledge, 
not tabulated. Nevertheless, we can make analytical progress by realizing that we only have 
to take into account the states with m = 0, since the interaction conserves parity. The radial 
part of these states is given by 



(B2) 




it 



ia^)/y2 lowers the magnetic quantum number m of the angular 



where the operator 

momentum by one. The operator (a J. + ia) y )/\/2 raises this quantum number by one. Here, 

the operators a\ = \Jmuji/{2Ti){xi — ifi/{muji)) are the usual harmonic-oscillator creation 
operators. The ground state is denoted by |0). The creation operators commute and hence 
we can rewrite the radial wave function of the state as 



n=0 



n\(n r — n)\2 r 



(B3) 



where |n)j denote the normalized eigenstates of the hamiltonian Hi = p 2 /(2m) + muJix]/2 
of the one-dimensional harmonic oscillator. In the derivation of this expression we used 
Newton's binomium to rewrite the n r -th powers of the operators in the right-hand side of 
Eq. (|B2|) as a sum of n r terms. With this result, the normalized wave functions of the excited 
states of interest are given by 



X 



n r ,n z 



E 

.n=0 



n\(n r 



n 



:e-^ + y 2 ^H 2n (x/l r )H 2(nr . n) (y/l r : 



xe- z2 /^H nz (z/h 



(B4) 



The overlap integrals of two excited states of this form with a gaussian of arbitrary width 
are tabulated EjJ. The final result for the matrix elements is given by 



Vn 



r ,n z ;m r ,m z 



(t) 



N c {t) 
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1 



n r m r min(2n,2m) min(2(n r — n),2{m r —m)) 

EE E E 

n=Om=0 £i=0 1=0 



1 



{2n)\ 2(n r -n) 



n\(n r — n)\ m\(m r — m)\ (2n — k)\ (2(n r — n) — /)! 
2m \ / 2{m r — m) 



x 



k 



I 



(2(m + n — k) - 1)!! 



x (2(n r + m r -n-m-l)- 1)!! ((q r (t)/l r ) 2 + 1)*+'-^-^ 



x 



min(n z ,m z ) 

^ ^!(n 2 -g)!(m z -g)! 



x((g z (t)// 2 ) 2 + iy/^ + n z )/4 {nz + mz _ 2q _ !),, 



(B5) 



if n z + m 2 is even, otherwise it is equal to zero. 
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Fourier transform of the gaussian ansatz, and assumed the energy of a condensate atoms to 
be negligibly small with respect to the energy of an ejected atom. Moreover, in the retarded 
propagator of the ejected atoms we neglected the mean-field energy with respect to the kinetic 
energy of the ejected atoms. Here we go beyond these approximations. 
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